Physical evolution in Loop Quantum Cosmology: The example of vacuum Bianchi I 
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We use the vacuum Bianchi I model as an example to investigate the concept of physical evolution 
in Loop Quantum Cosmology (LQC) in the absence of the massless scalar field which has been 
used so far in the literature as an internal time. In order to retrieve the system dynamics when 
no such a suitable clock field is present, we explore different constructions of families of unitarily 
related partial observables. These observables are parameterized, respectively, by: (i) one of the 
components of the densitized triad, and (ii) its conjugate momentum; each of them playing the 
role of an evolution parameter. Exploiting the properties of the considered example, we investigate 
in detail the domains of applicability of each construction. In both cases the observables possess 
a neat physical interpretation only in an approximate sense. However, whereas in case (i) such 
interpretation is reasonably accurate only for a portion of the evolution of the universe, in case (ii) 
it remains so during all the evolution (at least in the physically interesting cases) . The constructed 
families of observables are next used to describe the evolution of the Bianchi I universe. The 
performed analysis confirms the robustness of the bounces, also in absence of matter fields, as well 
as the preservation of the semiclassicality through them. The concept of evolution studied here and 
the presented construction of observables are applicable to a wide class of models in LQC, including 
quantizations of the Bianchi I model obtained with other prescriptions for the improved dynamics. 

PACS numbers: 04.60.Pp,04.60.Kz,98.80.Qc 



I. INTRODUCTION 

Loop Quantum Cosmology (LQC) [1] is an area of 
physics which deals with the quantization of symmetry 
reduced gravitational systems by adopting similar meth- 
ods to those employed in Loop Quantum Gravity Q. 

Not surprisingly, the kincmatical structure underlying 
LQC was rigorously established for the first time for the 
simplest of all the cosmological models, namely, those 
describing a homogeneous, isotropic, and spatially flat 
geometry 0). The analysis of the dynamics of such a 
kind of geometry in the presence of a massless and min- 
imally coupled scalar field has shown that the big bang 
singularity is resolved dynamically and is replaced by a 
quantum bounce 0,i]. LQC has been further extended, 
with diverse levels of rigor, to other similar models with 
different topology ((J or nonvanishing cosmological con- 
stant Q, and furthermore to more general settings such 
as anisotropic systems || H, E3], or even to inhomoge- 
neous situations [III • The robustness of the singularity 
resolution features has been confirmed within an exactly 
solvable version of LQC m E3, 

and the mathe- 
matical foundations of its elements have been discussed 
15, 16]. In addition to the studies of the genuine quan- 
tum theory, there exists an extensive amount of work 
at the level of the effective classical dynamics [13, fl8| . 
which provides important insights into the properties of 
the quantum geometry in cosmological scenarios 
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So far, all the analyses of the dynamics carried out in 
the genuine quantum theory have employed partial ob- 
servables parameterized by an emergent time. In most 
of the models considered in the literature 0, [E H S H> 
ITol . H3 , [3], the role of time was played by a massless 
scalar field, which (unlike the geometry degrees of free- 
dom) was quantized adopting a standard Schroedinger- 
like representation. However, in many cases (vacuum 
homogeneous universes, black hole interiors) such pos- 
sibility is not at hand. In this article, we overcome this 
difficulty by constructing and studying the properties of 
various families of unitarily related observables parame- 
terized by the geometry degrees of freedom. As an ap- 
propriate test bed for our analysis, we study a model of 
a vacuum homogeneous universe of the Bianchi I type. 

The Bianchi I system itself has been extensively ana- 
lyzed in the literature. After a preliminary analysis of 
the kinematics [2l| . an attempt to quantize the model 
was made in Ref. Q, and the corresponding effective 
dynamics was studied in Ref. [22| . There, however, a 
homogeneous massless scalar field was coupled as matter 
content and identified as an internal time. 

More recently, the system in vacuo has been thor- 
oughly quantized [9(. In that work, the studies were re- 
stricted to the model with the compact three-torus (T 3 ) 
topology, in order to provide a basis for further studies 
of a much more general system, namely the inhomogc- 
neous Gowdy T 3 model with linear polarization [Til ]. In 
Ref. @ both the physical Hilbert space and a complete 
algebra of observables were constructed; nonetheless the 
dynamics of the system was not investigated. 

Let us comment that the extension of a quantiza- 
tion prescription known as improved dynamics from the 
isotropic to the anisotropic setting in homogenous LQC 
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has given rise to the proposal of two concurrent prescrip- 
tions for Bianchi I. In this article, following Ref. Q, 
we use the prescription proposed in Ref. For non- 
compact models and just in the form presented in that 
paper, that prescription is not invariant under certain 
transformations [23(, but for compact models, including 
the one considered in Ref. 9], the prescription is free 
from this drawback. The desired invariance is recovered 
in noncompact situations with an alternative prescrip- 
tion, which was recently developed and motivated in Ref. 
(To| appealing to the hypothesized relationship between 
the degrees of freedom of LQC and full LQG. For our 
discussion, nonetheless, the following reasons explain our 
use of the former of these prescriptions, (a) The goal of 
our work is a methodological development of a formalism 
of unitary evolution, defining it in a precise manner with- 
out introducing an additional matter field as a clock. For 
that purpose, the model of Ref. @ is a better candidate, 
as the structure of the Hilbert space and the properties 
of the states are known in detail and have been presented 
in the literature, something which is necessary to estab- 
lish the link between physical intuition and the exact 
mathematical implementation of our method. (b) The 
specific properties of the model of Ref. [9j allow to study 
the limitations of applicability of some of the constructed 
methods (see in particular Fig. [3]and the related discus- 
sion in Subsec. IV E[) . limitations which would have been 
missed had one used the model of Ref. [ijj ■ (c) Our con- 
struction, although applied here to a concrete model, is 
intended to be reasonably universal. In particular, it is 
almost directly applicable to other possible prescriptions 
for the quantization of the Bianchi I model, including 
that of Ref. [l(| (see the discussion in Sec. IVIip . 

Employing the knowledge accumulated in previous 
works about the properties of the selected model, and the 
explicit and complete quantization chosen for it, we will 
analyze here two possible constructions of families of uni- 
tarily related observables, which are respectively param- 
eterized by: (i) a coefficient of the densitized triad, and 
(ii) the variable conjugate to it. For the sake of clarity in 
the presentation, as well as for comparing the predictions 
of geometrodynamics and LQC, we will apply the con- 
struction to both the LQC version of the model and its 
Wheeler-DeWitt (WDW) counterpart. We will see that, 
whereas in a WDW-type quantization our construction 
with respect to both of the specified parameters provides 
observables with a precise physical interpretation, in the 
loop quantization scheme this property is achieved only 
in an approximate sense. Furthermore, the accuracy of 
the interpretation depends only on properties of the state 
in case (ii), while for (i) the dependence is also on the 
evolution epoch. 

The paper is organized as follows. Sec. [II] is a review 
of the LQC quantization of Ref. @. In Sec. IIIII we 
perform the WDW quantization of the model and develop 
an interpretation for the evolution of observables for the 
two considered choices of internal time. In Sec. IIVI we 
show that the LQC states converge in a certain limit to 



a combination of WDW states. The evolution picture 
in the LQC theory is analyzed in Sec. [V] and Sec. [VI] 
for the choices of internal time (i) and (ii), respectively. 
In Sec. I VIII we discuss the main results of this work. 
Besides, two appendices are added. One of them extends 
the discussion on the WDW regime to take into account 
the union of different sectors. Finally, a brief description 
of the effective classical dynamics, associated with the 
genuine LQC dynamics, is given in Appendix [B] 



II. THE FRAMEWORK 

In this section we briefly review the basic features of 
the polymerically quantized vacuum Bianchi I model. It 
is a summary of the main results obtained in Ref. 0, 
which the reader can consult for further details. 



A. Classical theory 

The Bianchi I model represents spatially flat and ho- 
mogeneous spacetimes. We will consider the case of a 
compact topology: that of a three-torus. The spacetime 
metric can be written in the form 12211 
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(2.1) 



where {dx 1 } is the fiducial co-triad, N is the lapse func- 
tion, and pi/ (Air 2 ) are the nontrivial components of the 
densitized triad in a diagonal gauge. The correspond- 
ing components of the Ashtekar connection are c 1 /(2tt), 
such that {c l ,Pj} = SirGjSj. Here G is the Newton con- 
stant and 7 is the Immirzi parameter (that we assume 
positive for simplicity). Owing to homogeneity, the only 
constraint present in the model is the Hamiltonian one, 
given by [221 ] 

2 c 1 pic 2 p 2 + cVic 3 p 3 + c 2 p 2 c 3 p 3 
<-BI = o 77 — U ' ( Z - Z ) 
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V 



where V = y / \PiP2P3\ is the spacetime volume. Clas- 
sically, any of the triad coefficients jk is a monotonous 
function of the time coordinate t [22|. Moreover, the 
spacetime presents a curvature singularity at initial time 
t = 0, at which the universe stretches as an infinitely long 
line. 



B. Loop Quantum Cosmology kinematics 

In LQC the basic configuration variables are holono- 
mies of connections. The holonomy along an edge of ori- 
ented coordinate length lixpn in the direction i is defined 
as hf^c 1 ) = e MicVi , where 7j are the SU(2) generators 
proportional to the Pauli matrices, such that [t^Tj-] = 
e.ijkT k . The configuration algebra Cylg is the algebra of 
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almost periodic functions of c 1 , which is generated by 
the matrix elements of the holonomies N^d) = 



In the momentum representation, the states defined by 
these matrix elements are denoted by \fii). The com- 
pletion of the algebra Cyl s with respect to the discrete 



inner product 



for each fiducial direction 



Kin' 



provides the kinematical Hilbert space Jimn 
The elementary operators are the operators pi associated 
with fluxes, which are diagonal on the basis states \fii) 
of T^xin' and Ha 1 . , whose action shifts the label fa of the 
considered basis by 

In order to express Cbi in terms of holonomies, instead 
of connections components, one introduces the curvature 
tensor associated with the Ashtekar connection and de- 
fines it as in gauge lattice theories. With this aim, one 
considers a loop of holonomies and shrinks its area to its 
minimum value, which is nonzero owing to the discrete- 
ness of the area spectrum in LQG. To determine the mini- 
mum 1/Kfa for the coordinated length of the holonomy we 
assume (as in Ref. Q) square loops of fiducial area 4-7r 2 /i 2 
with minimum physical area equal to A = 2\/37rjlp l , 
with Zpi = V Gh being the Planck length. This implies 
that p%\pi\ — A in the sense of operators [HI- On the 
other hand, one appeals to the so called "Thiemann's 
trick" to define a well-defined inverse volume operator 

[1/V]. As a final result one obtains a quantum Hamil- 
tonian constraint Cbi expressed in terms of our basic 
operators and which is well-defined in the domain Cylg 

Si- 

Owing to the relation /uf \pi\ = A, the shift produced by 
the operator A/±p i on the basis states \ fj,i) is not constant. 
One can nevertheless relabel these states with an affine 
parameter 



Vi(fj,i) = 2=3 * 8gn.(fii)\m\* , 



(2.3) 



to make the shift uniform [H, 0| . The action of the basic 
operators in the relabeled states is Ai^Jvi) = \vi ± 1) 
and Pi\Vi) = S^Asgn^H 2 / 3 !^). 

We have symmetrized the resulting Hamiltonian con- 
straint Cbi in such a way that it annihilates the sub- 
space spanned by the "zero- volume" states, i.e. the states 
®i\vi} with any Vi equal to zero, and leaves invariant its 
orthogonal complement. For the nontrivial solutions to 
the constraint, we can thus restrict our study to that 
complement. Owing to this restriction the singularity is 
resolved at the kinematical level @. 

It is convenient to densitize the Hamiltonian constraint 
to find its solutions. The densitized Hamiltonian con- 
straint is given by the operator 
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(2.4) 



which is well-defined in the considered domain (linear 
span of tensor products of states such that none of 

the Uj's vanishes). The bijection = (ip\ [1/V] 1 ^ 2 in the 



dual of that domain relates the nontrivial solutions (ip\ 
of Cbi to the solutions (tp\ of its densitized version Cbi- 
This operator turns out to be given by 



C 



BI 



r 



eie 2 + 6163 + e 2 Oj 



(2.5) 



where the operator O^ is symmetric in the domain 
spanned by the basis states \vi) (wj ^ 0). Its action is: 



Oi\vi) 
where 



A 



i—[f + ( Vi )\ Vi + 2) - f-(vi)\vt - 2)], (2.6) 



f±{vi) = g{vi ± 2)s±(v t )g(vi), 
s±(Vi) = sgn(w 4 ± 2) + sgn(ui), 



9(vi) = 




if Vi^O, 
if v, = 0. 



(2.7) 
(2.8) 

(2.9) 



Let us note that the operator O^ leaves invariant the 
Hilbert subspaces Tif. defined as the Cauchy completions 
with respect to the discrete inner product of the spaces 

Cyl± - span{K); Vi € C%}, (2.10) 

where Cf. are the semilattices of step two 

C% ={±(^ + 270,71 = 0,1,2,...}, e» 6(0,2]. (2.11) 

Therefore the constraint Cbi superselects the kinemat- 
ical Hilbert space in different separable sectors. If we 
choose, for instance, a positive orientation for the triad 
in the three fiducial directions, we can restrict our study 
to the kinematical Hilbert space Jit = ®iJi£., with 

£*= [E\, £2, £3)- 

C. Physical Hilbert space 

In order to obtain the physical Hilbert space, it is 
enough to analyze the spectral properties of the operator 
Qi. It can be proven essentially self-adjoint with domain 
the dense set Cyl^ . In consequence, Cbi is essentially self- 
adjoint in the invariant domain Gy\t = <g>jCyl^. This 
property allows us to apply the group averaging method 
[l^lll] to determine the physical Hilbert space W? hy . 

In addition, the spectrum of O.; is absolutely contin- 
uous, coincides with the real line and is nondegenerate. 
Furthermore, the coefficients e e J. (2n + £,-) (n G N + ) of its 
generalized eigenfunctions with eigenvalue turn out to 
be determined just by the initial data efj.(£i) [see Ref. 
[9j for the explicit expression -Eq. (45)-] . 

These eigenfunctions are formed by two components, 
one of them with support in the semilattice of step four 

Wr+ ={(£i + 4n),n = 0,l,2,...} 
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and the other in the displaced semilattice 

{i) Cl +2 = {(ei + 2 + 4n),n = 0, 1, 2, ...}. 

Both of these components are generalized eigenfunctions 
with eigenvalue w 2 of the operator 0?, which is propor- 
tional to the gravitational part of the densitized con- 
straint in the isotropic case. These components have 
a relative phase of ±7r/2 and, hence, e £ j\vi) oscillates 
rapidly when Vi varies in the semilattice Cf. . 

Applying the group averaging procedure, we obtain the 
following form for the wave function of the physical states 



= J ^ dw 2 dw 3 $(^2, u 3 ) e^ l(w3)ta , 3) (ui) 

x e e J 2 (v 2 )ell(v 3 ), (2.12) 



with 



UJ 2 + &3 



(2.13) 



Here $(072,^3) belongs to the physical Hubert space, 
which turns out to be 



(2.14) 



In what follows, unless otherwise specified, uj\ will be 
the function U)\{w2,uz) given in Eq. (|2.13|) . 



III. WHEELER-DEWITT ANALOG 

In this section we will study the quantization of the 
vacuum Bianchi I model in the WDW approach formu- 
lated in terms of the connection coefficients, which will 
be treated in a standard (non-polymeric) way. 

A. Kinematics and scalar constraint 

Like in the loop quantization, we will work in the 
triad representation. As kinematical Hilbert space of 
the WDW quantization we take 2£icin = ®i2HKm w ^h 
21kin = L 2 (R, dvi) [13]. The measure is the usual 
Lebesgue measure, and not the discrete one. In this rep- 
resentation the operator pi acts by multiplication by the 
factor pi = 3 1 / 3 Asgn(wi)|«i| 2 / 3 , just as in the loop quan- 
tization. On the other hand, we promote the connection 
coefficients to derivative operators, 



& = &/*2\v i \ 1 ' 3 d Vl , 



(3.1) 



to preserve the Dirac rule [c 1 ,^] = ih{c l ,pj}. 

We denote by O^ the counterpart of the classical quan- 
tity c l pi in the WDW quantum theory [defined on the 
Schwartz space <S(R)]. We choose for 0^ the symmetric 
factor ordering which is analog to that used in the loop 
quantization, to simplify the comparison with it: 



4 = iVZA^/\iH\[sgi\{vi)d v 



(3.2) 



It is well defined in the distributional sense and can be 
rewritten in the simpler form 0^ = iy3A(l + 2vid Vi ), 
where we have disregarded the noncontributing term 
|i?i|<5(Ui). Then, the WDW quantum counterpart of the 
classical densitized Hamiltonian constraint can be con- 
structed exclusively from the operator O^ in the following 



- BI _ 



©102 + 0103 + 0203 



(3.3) 



The operator 4 is essentially self-adjoint on 2ikin- 
Furthermore, its restrictions to each of the subspaces 
2ikin = i 2 (R ± ,c?w i ) are essentially self-adjoint as well. 
On each of these subspaces, the spectrum of 4 is ab- 
solutely continuous, coincides with the real line and is 
nondegenerate, as happens to be the case with its analog 
Qi in the loop quantization. Moreover, its generalized 
eigenfunctions, with generalized eigenvalue u>i, have the 
form 



1 



y / 2na\v i 



■. exp 



In \v. 



(3.4) 



127T7^pj. They provide an orthonor- 



where a = 2V3A 

mal basis for 21kin- ^ n ana logy with the procedure fol- 
lowed in the loop quantization, we will restrict the study 



to U£ in 



iU 



';+ 

Kin' 



B. Physical Hilbert space and observables 

In order to obtain the physical Hilbert space, we apply 
the group averaging method like in the loop quantization. 
Obviously, the physical Hilbert space obtained is 



W Phy 



\uj 2 + 0J3\daj 2 duJ 3 



(3.5) 



The wave function of the physical states has the following 
form: 



/ duj 2 duj 3 #(w 2 , ^3) e^iCwi) 
Jr 2 



(3.6) 



with $(w 2 ,^3) e 

As in the LQC model @ , the operators O a (here a = 
2,3), which multiply the wave function by u> a , together 



with 



-l\U> 2 + W 3 | 



-l/2fl 



1^2+^3 



1/2 



provide a complete 



set of observables, which are essentially self- adjoint op- 
erators on the domain 6>(R 2 ) C 7i Phy . Nonetheless, 
while the first pair are Dirac observables which corre- 
spond to classical constants of motion, the second one 
are not. As a consequence, this set turns out not to be 
adequate to introduce a nontrivial concept of evolution. 
In the next subsection, we will construct another pair of 
Dirac observables that will complete the set formed by 
a (a = 2, 3), and in terms of which we will be able to 
develop an interpretation for the notion of evolution. 
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C. The evolution 

In any gravitational system, as the one considered 
here, the Hamiltonian is constrained to vanish. There- 
fore, there is no well defined notion of evolution in the 
model. Nonetheless, one can try to select, in the con- 
figuration space, a coordinate t (an internal time) and 
define a map between Hilbert spaces 2£ t , where t takes 
values in a certain set U and H t is the space of "initial 
data" given by the restriction of the wave function to 
the surface t = const. If there exists a unitary transfor- 



mation P t : n Phy 



Tt t , then each i-slice contains all 



the information needed to determine the physical state 
(i.e. the system is closed). Furthermore, if, in addition, 
the "identity" map between the spaces H t (given by the 
trivial identification of data at different times) is also uni- 
tary, one can define a unitary evolution in 7£ t . This is 
achieved by composing the inverse transformation P^ 1 , 
a transformation in the family P_ t (for a different value 
of t) , and the above identification of data. 

In our case, since classically any of the three triad com- 
ponents is monotonous along all the dynamical trajecto- 
ries, one can select one of the Vi's as an internal time. 
Taking into account that in the description of the phys- 
ical states p.6p we have already eliminated uj\ in terms 
of u 2 and oj 3 , it is then most natural to choose t := v\. 

We can easily introduce the initial data spaces labeled 
by v\ = const (which we will also call "slice" spaces) via 
the transformation 



® Vl (V2,V3) = / dw2dW3$ (w2,W3)e^ 2 (^2)^3(^3) 
JR 2 

■=^2,(^3)^^1), 

with $ (u>2, 0J3) belonging to 

H V1 = L 2 (R 2 , \u 2 + W3 ||^>i)r 2 d W 2d^ 3 ) 
2 2 

= L (R ,2iravi\uj 2 + uj 3 \duj 2 dio 3 ). 

Since given in Eq. (|3.4[) . never vanishes, all the 

slice spaces H Vl are unitarily related to 7£ Phy ■ 

Introducing an additional rescaling one can define the 
alternate transformation 



§L Vl (V2,V 3 ) 



y/2iravi 
§^(^2,^3) := P^<|(w2, ^3) 



dw 2 dw 3 , \ t \ r \ 



:= $(w 2 , t^v^Travie^ (wi) G 2£ Vl > 



where 



pi . -ryPhy 'yji 
—v\ • LL J-kvi 

r 2/TD)2 



K V1 = L 2 (R',\LU 2 +LU 3 \du 2 dL0 3 ). 



(3.10a) 
(3.10b) 



Since P[ Ul is unitary and the spaces H Vl coincide (so that 
the corresponding map identifying states at different val- 
ues of v± is unitary as well), it is clear that P v leads to 
a unitary evolution for the system. 



A more elaborated manner of providing a concept of 
evolution is by introducing a family of partial observables 
[HI, [29| , which are to be related via unitary transforma- 
tions. A natural construction of partial observables from 
the kinematical ones is available if there exists an internal 
time which provides a unitary map between initial data 
spaces. On the formal level, one can do this through the 
group averaging procedure. Alternatively, if there ex- 
ists a decomposition HKin —l£®2£ (where is some 
Hilbert space of functions depending on t only)[3(| and 
there is an operator O 1 : H' — > H' , one can build Ot 
(which measures the corresponding quantity "at a given 
time t") as an operator whose action is defined through 
the following sequence of operations: 

(i) Take an initial data slice ipt G Ht corresponding to 
the value t of an emergent time. 

In many cases the spaces H t differ from 7{ . There- 
fore one has to define a transformation 7i t — > H! 

a _ _ 

(ii) Then act with the "kinematical" observable O' cor- 
responding to the measured quantity, and 

(iii) use the relation between physical and slice Hilbert 
spaces to find the element of 71 Phy corresponding to 
the result of $u§. 

In the model under study, the above method can be 
applied for example to the kinematical observable ln(w a ) 
(where a = 2,3), which acts on elements of ^Kin as a 
multiplication operator. As the internal time coordinate 
we still choose vi, so that H t '■= tL v ; given by Eq. (|3.8|) . 
The space H' is the product 



L 



-) 2 ,dv 2 dv 3 ) (3.11) 



Since 7i v and 71' are different, we introduce a uni- 
tary transformation 2£„ — ► 2£- In the corresponding 
w a -representations, it is given by the map 



in 



:= V27tcwi|cj 2 + W3| 2 j^j (lu 2 ,uj 3 ). 

In the ^-representation, this transformation is simply 
the map $„ 1 (v 2 ,v 3 ) -> x. Vi (v 2 ,v 3 ), with $„ 1 (v 2 ,v 3 ) given 
in the first line of Eq. (|3.7|) and 



X Vl (v2,v s ) 



duj 2 du 3 x vi (w 2 , ^3)e W2 (^2)^3 («3)- 



(3-13) 

Employing the transformations between the intro- 
duced Hilbert spaces, we finally obtain a family of ob- 
servables hi(v a ) Vl acting on the physical Hilbert space, 
interpretable as "the value of ln(i> a ) at the fixed time 
ui". The operators are well defined on the Schwartz 
space 5(R 2 ). Their action in this domain, which can be 
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deduced taking into account that they act by multiplica- 
tion on the corresponding kinematical states x i v 2, u 3 ), 
turns out to be 



[ln(u a )„ 1 $](cJ2,w 3 ) = 
x d. 



UJ 2 +UJ 3 \ 



LU 2 + w 3 |?$(w 2 , uJzj^ivi) 



The observables ln(w ) Ul , together with the constants 
of motion OJ^ := © a , form a complete set of Dirac 
observables. Whereas within each family (corresponding 
to a = 2,3 respectively) the observables a | w do not 
change with v\, the observables ln({i a ) 1)1 do not coincide, 
and are related at different times v± and v* via an oper- 
ator Q ... : H Phy -» K Phy such that 



,|>](w 2 , w 3 ) 



-!— — J>(w 2 ,w 3 ■ 3.15 



uT e 



The form of the eigenfunction (13.41) implies immediately, 
that these operators are both invertible (Q . = Q . ) 

and unitary on 7£ Phy . Therefore, the relation between 
observables at different times, 



ln(# )„* 



(3.16) 



is unitary. 

As a consequence, the families of observables ln(z) a ) Ul 
define on the physical Hilbert space 7Y Phy a unitary evo- 
lution that is local in the emergent time v\. In contrast, 
as we will see in Subsec. IV B| the direct application of 
the above construction in the loop quantization does not 
lead to a unitary evolution since, in that case, the analogs 
°f Q v v * fail t° De unitary operators. Nonetheless (as we 
will see in Subsec. IVI Bj) families of unitarily related ob- 
servables can be defined once we use, instead of v\, its 
conjugate momentum, denoted by b\, which provides a 
suitable emergent time in the loop quantization. 

To compare the dynamics predicted by the constructed 
families of observables with the classical dynamics, let us 
calculate the expectation values on some class of states 
which are semiclassical at late times, namely Gaussian 
states peaked around large values of lo\ and u 3 : 



J>(w 2 , LU 3 ) 



K 



y/\u 2 +U) 3 \ a 



n 



71 e i0 a u„ 



, (3.17) 



where the constant coefficients A a and B a are 

A a = \\uj 1 {Lu 2 ,Lu 3 )Lu- 1 <t\\ 2 , (3.19a) 
B a = a(§l\u2 + u%\-*{-id Ua )\w2+ua\*&). (3.19b) 



Evaluating them for the Gaussian form (13.171) of $ and 
taking the limit a a — > (for both a = 2,3), we obtain 
the following trajectory 



(£IM*a)„i£) = 



lni)i+a/3 a (3.20) 



which agrees with the classical one. This result implies 
that in the WDW theory the singularities of the vacuum 
Bianchi I universe are not resolved dynamically (in the 
sense of the trajectories defined by the expectation values 
on semiclassical states). Actually, the lack of singularity 
resolution is a general property of all the states for which 
the coefficients A a and B a defined in Eq. (|3.19[) arc finite. 

To analyze the behavior of the dispersions (A ln(v a ) Vl ) , 
we first find the expectation values of lii 2 (v a ) Vl in a way 
similar to the derivation of Eq. (|3.18p . They read 



(tlln 2 ^)^) = W„ In 2 vi + F a hi + X a , (3.21) 



where 



W a = \\Lol(Lo 2 ,u; 3 )LU- 2 §>\\ 2 , (3.22a) 
Y a = -2ia($||w 2 +w 3 | - *wi(w2,W3)a;~ 1 '(c?wJ 

X |w 2 + Uz\*uJi{u2,uz)uJa l &l (3.22b) 
X a = -a 2 (l||w 2 + w 3 |-5a^ |w 2 +w 3 |5j>). (3.22c) 



In Eq. (|3.22b[> . there is no summation over the in- 
dices a. Using the standard relation (Aln(-0 Q )„ 1 ) 2 = 
(ln^n)^ ) — (hn.(va) Vl } 2 we can easily find the dispersions. 
In particular, we see immediately that once these disper- 
sions are finite in some epoch, they remain so throughout 
all the evolution. Furthermore, for states for which the 
expectation values B a and X a [defined in Eqs. (|3.19p 
and (|3.22[) ] are finite, the relative dispersions approach 
constant values in the large v\ limit, values which are 
determined by the relative dispersions of w 2 (cij 2 , ui 3 )to~ 2 : 

Hm (A\n(v a ) Vl ) _ (A[u; 2 (o; 2 ,u; 3 )^- 2 ]) 



vi^oo (ln(Va) Vl ) 



(u((u;2,UJ 3 )u)a ' 



Here we have used the shorthand (O) := ($|0 $) for any 
operator O. 



where K is a normalization factor such that II $11 = 1, 



and the factor \u)2 + w 3 | 



-1/2 



compensates the nontrivial 



factor in the measure of the physical Hilbert space (|3 . 5|) . 

For a general state $, using directly the explicit form of 
the observables (I3.14| and integrating the inner product, 
we find that 



(*| hi^a)^!) = A a \n Vl +B a 



(3.18) 



D. 6-representation 

In the loop quantization, the evolution with respect to 
the internal time v\ , analog to the one constructed above 
for the WDW theory, fails to be unitary. Nevertheless, 
the use of the momentum conjugate to v\ as internal 
time provides a good notion of unitary evolution. Let us 
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study this choice of time in the WDW quantization as 
well, both for completeness and in order to introduce the 
procedure in a simple setting, where the difficulties inher- 
ent to the more complicated nature of the LQC model are 
absent. 

Given the classical variable introduced in Eq. (12. 3|) . 
one can define the conjugate momentum 

h := VAc l \ Pl \~K (3.24) 

such that 

{b i ,v J } = 25 ij . (3.25) 

In the WDW quantum theory, the unitary transforma- 
tion between the "position" and "momentum" represen- 
tations is given by the Fourier transform: 

t£VK&i) = 1T7= I dv^iv^e-^. (3.26) 

It is worth emphasizing that this is a unitary transforma- 
tion from the Hilbert space of square integrable functions 
in the ^-representation to the same Hilbert space in the 
^-representation: 

£ ■ 2fc in = L 2 (R, d«0 -» fC mn = L 2 (R, dh) (3.27) 

Under this transformation the elementary kinematical 
operators transform as 

Vi -> 2id bi , d Vi -> ibi/2, (|3~28|) 

where bi acts in the new representation as a multiplica- 
tion operator. Hence, the transformed of 0^ is 

= -i\/3A(l + 2bid bi ), (3.29) 

which coincides with the original operator 0^ (up to 
a sign), both being defined in identical Hilbert spaces. 
Therefore, in the WDW quantum theory, working in the 
^-representation is completely equivalent to working in 
the Wi-representation. In particular, we can regard b\ 
as the internal time, change the representation only in 
the direction 1 , and define the Hamiltonian constraint in 
the kinematical Hilbert space fj^ltn ®a i£(L ( a = 2,3), 
where = L 2 (R+,dbi) [13], by replacing the opera- 

tor ©! with the operator £(©i) in Eq. (|3.3p . Hence, we 
can repeat exactly the construction introduced in Subsec. 
|HICJ substituting v\ by b\ and e^Wi) by e_ Ul (bi). 

Whereas the ^-representation does not introduce any 
novelty or advantage in the WDW quantization in com- 
parison with the ^-representation, we will see in Sec. IVII 
that there is a big difference between both approaches in 
the LQC quantization. 

IV. WHEELER-DEWITT LIMIT OF THE LOOP 
STATES 

The comparison of Eqs. (|3.5j) and (|3.6|) with Eqs. 
(I2.12p - (|2.14[) shows that the physical Hilbert spaces of 



the LQC and WDW quantizations, as well as the struc- 
ture of the corresponding wave functions, are identical. 
The difference between both quantizations is captured in 
the different form that the eigenfunctions of the opera- 
tors Oi and 0j possess. On the other hand (as it will 
be shown in Subsec. IIV A| . the LQC eigenfunctions con- 
verge for large Vi to some combinations of their WDW 
analogs. This feature allows one to regard the dynamics 
of the LQC universe as certain form of "scattering" of 
WDW states, incoming from a distant past and outgoing 
to a distant future [H[ . Mathematically, this behavior is 
described by the analog of a scattering matrix p s acting 
on the incoming WDW state: 

|J>)out = p s |l)in, (4.1a) 

x ,02(^2, <4)p3(a;3,c4) (4.1b) 
pi(u)i,u}-) := (e^Jpile^). (4.1c) 

In turn, the matrices pi are determined by the WDW 
limit of the eigenfunctions of 0^. The exact form of this 
limit will be investigated in Subsec. IIV Al The result 
will be applied in Subsec. IIV Bl to describe the limit of 
physical states. Finally, Subsec. IIV CI deals with the 
effect of the commented scattering on the dispersion of 
observables. 



A. Limit of the eigenfunctions 

Let us restrict ourselves to one superselection sector 
(defined in Subsec. Ill C|) . e.g. that corresponding to 
functions supported on the lattice £+. The eigenfunc- 
tions efj. of 0^ with eigenvalue u>i were provided already 
in Ref. [§] [see Eq. (45)], although in the form presented 
there one cannot easily determine their large Vi behav- 
ior. In order to find it, we have to analyze the operator 
©i itself. Its properties were discussed in detail in Ref. 
Q as part of the proof of its self-adjointness. First, all 
of its eigenspaces are one-dimensional, and the complete 
solution is determined by an initial value efj.(ei). In par- 
ticular, the freedom in the choice of a global phase for the 
eigenfunctions is removed by demanding the positivity of 
this initial value e^j (£j). Second, as discussed in Subsec. 
Ill CI one can split the support into two subsemilattices 
( 4 >£f., with £j £ {£;,£.; + 2}. The restriction of efj. to 
each subsemilattice is an eigcnfunction of the operator 
0? with eigenvalue uif : 

Qlei\=^e% e% :=e^| (4)£ +. (4.2) 

The operator 0?, in turn, is a second order difference 
operator with real coefficients (similar in structure to the 
evolution operator defined in Ref. Q). It also has the 
property that all its eigenspaces are one-dimensional and 
the eigenfunctions are determined just by their initial 
value (at «j = e,). 
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To check the existence of the WDW limit of the eigen- 
functions, we implement the method used in Ref. [34j : 

(i) First, we represent the values of the eigenfunction 
at two consecutive points of ^Cf. by vectors i[)(v) 
of coefficients of its decomposition in terms of the 
WDW eigenfunctions e±i W4 i evaluated at this pair 
of points. 

(ii) Next, we rewrite the eigenfunction equation as the 
first order one acting on the vectors ip(v). It has 
a form of a 2 x 2 matrix [denoted from now on as 
B(v)}. 

(iii) Finally, we calculate the asymptotic expansion of 
the matrix B in the large v limit. 

An explicit calculation shows that the considered matrix 
is of the form B(v) = I + 0{v~ 2 ). Thus, there exists a 
well defined limit 

ip := lim tp(v). (4-3) 

v — >oo 

This immediately implies the convergence to a combina- 
tion of functions in the WDW basis, represented by the 
coefficient vector rp. The form of B implies also that the 
rate of convergence of ip(v) is at least of order 1/v. 

Actually, the fact that the operator 0| , after a suitable 
change of representation, differs from the isotropic evo- 
lution operator of Ref. [5| just by a compact term (see 
Ref. 0) allows us to apply here the numerical results of 
that reference. They show that the convergence is even 
faster, namely 

tjj(v) = ijj + 0(v-i). (4.4) 

Furthermore, from the reality of Of, it follows that the 
incoming and outgoing WDW plane waves contribute 
equally to the limit, that is 




To determine the normalization factor r we note that 
\im v ^ oc (d v e LJ .(v))/e LJ .(v) = 0. This fact, together with 
the sufficiently fast rate of convergence (|4.4|) . imply that 
the (kinematical) norms of both e E J. and its WDW limit 
(denoted here as e£*.) satisfy the relation 8||e^. || 2 = 
||e^. || 2 . As a consequence 

»•(<*) = V2z h (4.6) 

where Zi is a global phase which equals 1 for < 2 and 
— isgn(o;i) otherwise. 

The phase shift 4>(uJi) has a nontrivial dependence on 
u>i and needs to be found numerically. Luckily, the sim- 
ilarity of the operator O 2 with the evolution operator of 
an isotropic universe allows us again to apply directly the 
methods of Ref. [33| . The result is the following: 

<j)(uji) = (In |^| +a)(\ui\ +b) + c Si + R gl (\uj t \), (4.7) 



where a, b and cg i are constants, lim^^oo i^do;*!) = 0, 
and the dependence on enters only in the constant 
term and in the remnant part. Hence, for large |w;| the 
terms that affect the position and dispersion of the wave 
packet do not depend on the value of 

At this point, we can already write the exact form of 
the scattering matrix pi defined for just one subsemilat- 
tice. It reads: 

p i (u H ,Lu' i )=e- 2i 't>^6(u H +u;! i ). (4.8) 

Let us remember that the parts supported on differ- 
ent subsemilattices (individually of constant phase) are 
shifted in phase by tt/2. As a consequence, the common 
WDW limit for both of them does not exist (see Fig.[T]). 
Therefore, one cannot write an explicit form for the scat- 
tering matrix pi on the entire lattice Cj. . Nonetheless, 
the differences between subsemilattices manifest them- 
selves only through the constant phase shift and the 
remnant decaying for large As a consequence, when 
we consider the properties of asymptotic wave packets, 
peaked around large u> a , we can safely restrict the stud- 
ies just to one subsemilattice. 




FIG. 1: An example of eigenfunction of the operator Oi, 
corresponding to the eigenvalue u>i = 100 and the superse- 
lection sector gj = 2. The blue line (located on the imagi- 
nary plane) shows the part supported on the subsemilattice 
' 4 ^^. =4 , whereas the red line (real plane) is the part sup- 
ported on ' 4 )£g. =2 . 



B. Physical states 

The WDW limit of the eigenfunctions found in the 
previous subsection can be now applied in the analysis of 
the physical states. We start with a general state $, again 
restricting the study just to particular subsemilattices 
^Cf.- In order to find the limit of the corresponding 
wave function, we simply replace the basis functions efj. 
in Eq. (|2.12p with their large u, limits 

4 - r[e^ + e"^ e_J. (4.9) 
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Upon this replacement, the wave function is trans- 
formed into the function 

S2 , S3 =±i (4.10) 
x ^(0.2,0-3) ^2(^2) e^aC^), 
where s := (s 2 ,s 3 ), u} S {uj 2 ,uj 3 ) := ^1(32^2,53^3), and 



Sl=±l 



$g(u 2 ,u 3 ) = 2V2 J2 ^(5152^2, sis 3 w 3 )e 4Sl ^i ?) 

(4.11) 



,2_,2_,2 3 
z l z 2 z 3 



XS 



a=2 



where zf equals 1 for e i < 2 and —1 otherwise. Note 
that, in order to obtain the above expression, we have 
taken into account that u>i(—u>2, —U3) — — Wi(a;2, ^3) and 



r(s 3 uji) = s 



(l-*i)/2 



v22j (with u>i = u)g). 



If we compare Eq. (|4.10p with Eq. (|3.6p we see that 
each term, corresponding to a particular set of values 
of s a , has a form very similar to that of a WDW state 
with spectral profile $ = &g, the only difference being 
the replacement of the function with u>g in the in- 
dex of the basis functions. As a consequence, each of 
those terms can be (independently) considered as a state 
defined within a certain analog of the WDW quantum 
theory described in Sec. IIIII Such an analog inherits all 
the properties and structure of the original theory, except 
for the transformation of u>\ — ► u>g pointed out above. In 
particular, the inner product and the definitions of all 
the observables remain unmodified. 

This correspondence allows us to apply directly the 
definitions (|3.14[) of the observables \n(v a ) Vl to the terms 
(frg. Their expectation values (calculated on each term 
independently) satisfy an analog of the equation (|3.18|) 
with the direction coefficient A a replaced with 



A a ,g ■■= S a ||w s -(w2,W 3 )w a 



(4.12) 



To select the terms &g which actually contribute to the 
investigated limit we note that, if the state is localized 
(in the sense that it remains peaked around the trajec- 
tories defined by the expectation values of the consid- 
ered observables), only those terms corresponding to an 
s for which A a< g is strictly positive will have a signifi- 
cant contribution in the regime when the three v^s are 
all large. This requirement is satisfied only by the term 
with S2 = S3 = 1. 

The surviving term encodes the state of the genuine 
(untransformed) WDW theory. Therefore, the large vi 
limit of a localized state $ is simply given by a WDW 
state of spectral profile 

_ 3-z 2 -z 2 -~ 2 

$(^2,^3) = 2V2 ^2 s ~$(su 2 , sw 3 ) 

s=±l 

3 

xY[z i e u rt 3U *\ (4.13) 

i=i 



where u\ is again related to ui 2 and uj 3 via Eq. (|2. 13|) . 

Let us recall at this stage that the variable v\ plays 
the role of time; therefore it is proper to introduce a 
decomposition of the state into positive and negative fre- 
quency parts, corresponding to uj± > and u)\ < respec- 
tively. Physically, these components can be interpreted 
as moving forward and backward in time. An analog cor- 
respondence can be applied to other directions, defining 
the splitting into expanding and contracting components. 
Then, the change of sign in uj a corresponds to a parity 
reflection in x a = ln(w a ). As a consequence, the trans- 
formation {oJi} — > {— u>i}, which we have recognized as 
a symmetry of Eq. (|2. 13[) in our previous discussion, is 
the analog of the full PT (parity/time inversion) trans- 
formation. 

Let us now consider a LQC state with profile $(o;2, W3) 
which is localized in the sense explained above, and 
restrict our considerations to the corresponding WDW 
state with profile $(w2,w 3 ). For convenience, in the fol- 
lowing we will call 2v / 2$(w2, w 3 ) the reference state. The 
expectation values of \n(v a ) vi follow then the trajectory 
(j3.18p for some constants A a and B a . Furthermore, from 
Eq. (|4.13p , one can immediately see that the WDW state 
consists of two parts: one with the same parity/time ori- 
entation as the reference state (s = 1) and another which 
is PT reflected (s = —1). We denote them, respectively, 
by $ + and j>_ . If the reference state has a definite time 
orientation (only one sign of u)\ contributes) then the dis- 
tinguished components are, respectively, the part comov- 
ing in time with the reference state and the time-reflected 
part. Since the norms of both of these components are 
equal, any wave packet moving backward in time is fully 
reflected into a packet moving forward. This shows the 
presence of a bounce also in the internal time direction. 

In order to compare the trajectories of the expectation 
values of \n(v a ) in the considered components with those 
of the reference state, we note that, since their spectral 
profiles are related just by a rotation, and by the reflec- 
tion in uj a together with a possible change of sign in the 
case of the directional coefficients A a are the same 
for these component states and the reference state. The 
only change is in the parameter B a . As a consequence, 
the trajectories are just shifted with respect to the tra- 
jectory of the reference state. 

In the particular case of Gaussian reference states 
which are sharply peaked around large to*, 



$(^2,^3) 



K 



yj\u 2 +W 3 \ a 



n 



(4.14) 



(where K is chosen so that $ has unit norm), we can 
provide a quantitative estimation of the trajectory shift. 
Namely, for such states the phase rotation 4>(u>i) is well 
approximated by its first order expansion (with coeffi- 
cients D and E which depend on uj*) 



DK)=sgnK)(l + a + mK|), 



(4.15) 
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where a is the constant given in Eq. (|4.7p . In view of this 
approximation, one can almost directly repeat the calcu- 
lations of Subsec. IIII CI to find the formula equivalent to 
Eq. (|3.18[) . The resulting trajectories are 



pjav ](ui) = 



wi(w£,w£) 



[lnui -aD(uj\)) 



(4.16) 



•a[£>(wS)±/n 



Generically, these trajectories are disjoint, both between 
them and with respect to the trajectory of the reference 
state. Nonetheless, note that the two considered trajec- 
tories coincide in the case (3 a = 0. 



C. The dispersions 

In the discussion above, we have taken the limit a a — > 
(a = 2, 3) and considered the linear approximation to the 
variation of the rotation phase 4>(u>i). This essentially 
removes all the information about the behavior of the 
dispersions. However, in the scenario described in the 
previous subsection, where the LQC dynamics can be 
(in an asymptotic sense) viewed as the reflection of the 
WDW wave packet which moves backward in time 
and is contracting, into the moving forward and expand- 
ing wave packet $ , , it is important to ask how much 
the dispersion of J> + grows in comparison to the one of 

J>_ (or vice versa). For the isotropic model with mass- 
less scalar field, restrictive bounds on possible dispersion- 
growth have been found. This result can be extended by 
employing the exact triangle inequalities 33] which in- 
volve the dispersion in In \u> a \ and ln(t) a ) (for brevity, we 
suppress the subindex V\ in the latter of these operators). 

In the model considered here, given that the analysis 
of the WDW analog shows that the relative dispersions 
in \n(v a ) are the ones which approach constant (nonzero) 
values for large v\ [see Eq. (I3.23[) ]. we are interested in 
finding a weaker relation, involving exactly these relative 
quantities. 

We begin by recalling that the complete wave function 
<& is supported on the product of three semilattices , 
each of which is in turn the union of two semilattices 
of step four, ^C^^ (where ej <G {e^Ei + 2}). We then 
divide the support into eight sectors, corresponding to 
the eight products of the chosen semilattices of step four, 
and consider the respective restrictions of $ to each of 
them. In each sector, we apply the scattering scheme de- 
fined at the beginning of this section. In that scheme, the 
wave packet is transformed into $ + via the unitary 
rotation 



(4.17) 



i=l 



the reflection in the signs of u a , and a possible change of 
global sign. 



It seems reasonable to restrict our considerations to 
states such that the each of the components J>_ corre- 
sponding to the sectors introduced above is such that 
the coefficients A a , B a , W a , Y ai and X a , defined in 
Eqs. (|3.19[) and (|3 . 22[) . are finite. This ensures that each 
of these components has a well defined associated trajec- 
tory (|3.18p and a finite dispersion for every finite v\ . In 
addition, Eq. (|3.23p is satisfied. 

Neither the rotation (|4.17p . nor the reflection in the 
signs of Lu a , nor the possible change in sign of the wave 
function will change the expectation value or the disper- 
sion of the operator uif(ui2, ^3)a;~ 2 . Therefore, if the coef- 
ficients (|3~T9)) and (|33?2"1) are finite also for $ + , Eq. (|3~23|) 
implies immediately that the following holds 



(Aln(0) a ) 



(Aln(«) a ) 



(ln(«) a > + (\n(v) a ). 



(4.18) 



where ( • )± denote expectation values on $ ± , respec- 
tively. 

To prove that the finiteness of A a , B a , W a , Y a , and X a 
for J>_ implies the finiteness of these coefficients for $ + , 
we apply similar methods to those proposed by Kaminski 
and Pawlowski [33|. Since A a and W a are equal for both 
components, the only ones that require detailed analysis 
are B a , Y a , and X a . We first recall that 



B a 
X„ 



B 2 = 



(Av a y 



where 



V a = -ia\uj 2 + uj 3 \ 2 (d Ua )\dJ 2 +^3| 2 , 

V'a = 2W! (w 2 , W3)w a _1 P„ Wi (U) 2 , W3)w„ _1 ■ 



(4.19a) 
(4.19b) 

(4.20a) 
(4.20b) 



Knowing the relation between $_ and $ + , we can find 
the relations between their corresponding expectation 
values (|4.19p . They are 



3 

i=l 



(V a )+ = -<£>„)_ -2a} ^[duJi-Ui)])- (4.21a) 

-Ui)])- (4.21b) 



{V a ) + = -{V a ).-4aT i ('4[d t 



W a + < W a - + 2a^2{A[d Ua (f)(-uJi)])- (4.21c) 



i=i 

Actually, it is possible to estimate the terms related with 
dui a <p(—0Ji) employing that the function <p(ui) possesses 
the following properties [331 ]: 

\d u <ft{u)\ < Ci|mM|+Co, (4.22a) 
\uj8U{uj)\ < C 2 , (4.22b) 

where Co, C\, and C 2 are (positive) finite constants, 
which however may depend on the value of the subsemi- 
lattice label e$ and, in particular, may not have a global 
bound (in the whole interval of variation of this label). 
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These inequalities can be next used to relate the terms 
in Eq. (|4.2ip with the expectation values and dispersions 
of the operators ln|wi|. In order to do so, let us first 
define the multiplicative operators 



w a := <Ji(u>2, w 3 )w Q \ 
fiW :=<ln| Wl (a> 2 , W 3)|, 

4 B) -=<^al 



(4.23a) 
(4.23b) 
(4.23c) 



In the case of relation (|4.21a[) , the last term is bounded 
as follows 



(I4~24l) 



<Ci 



(|lnK||)_ + (|f7i 2 )|)_j+C (l + (^)), 
whereas for the term in Eq. (|4.21bp one has 

£(^>(-^)]> 

(|Si 2) |)- + (|^ 4) |>- +C [(wl 



i=l 



(I4~25l) 



<Ci 

Similarly, the sum in (|4.21cp satisfies 



£<A[^>(-^)]>- 



i=l 



(4.26) 



< C 2 



(Aln|w |)_ + <Anl- 



(2)\ 



Suppose now that in the state $_ the dispersions and 
expectation values of ln|o; Q |, and are finite, as 

(2) 

well as the expectation values of £a . This immedi- 
ately implies the finiteness of the right-hand side of Eqs. 
(|4~24l) - (|4~26|) (including that of (|fii 4) |)_). From this and 
the fact that Co , C\ , C2 < 00 it is straightforward to 
check that, if the coefficients B a , X a , and Y a correspond- 
ing to are finite, so are the ones corresponding to $ , . 
Therefore, we conclude that Eq. (|4. 18|) is indeed satisfied, 
as we wanted to prove. 

It is worth noticing that, in the previous discussion, 
the roles of $ + and can be interchanged, so that one 
can instead impose mild conditions of the type explained 
above on <& + and ensure then a good behavior for the 
relative dispersions corresponding to $_ . 

Finally, in Appendix [A] we show that our result (|4.18|) 
about the relative dispersions can actually be extended 
to the case in which one takes into consideration not 
just one isolated sector, but the whole ensemble of the 
eight sectors in which the LQC physical states admit 
a WDW limit. In this way, we arrive at the following 
conclusion. Consider a physical state described by the 
wave function <I> supported on the product of semilattices 
£+ x £+ x £+ . Suppose that the restriction of $ to the 
product of subsemilattices x ^Cf 2 x pos- 

sesses a WDW limit such that the component $ ± moving 



forward/backward in time has finite expectation values 
and dispersions for the operators 



ln(v a ) vi ,\n\aj a \,n^ ,wl. 



(4.27) 



Suppose also that this state has finite expectation values 
for the operators 



(4.28) 



Then 



(i) the corresponding component $ T moving back- 
ward/forward in time has also finite expectation 
values and dispersions with respect to the opera- 
tors (|4.27p . as well as finite expectation values for 
the operators (|4.28j) . and 



(ii) relation (|4.18|) holds for the ensemble of the WDW 
limits corresponding to all of the eight sectors de- 
fined by the restrictions to the different subsemi- 
lattices, constructed to reflect the relevant features 
of a complete LQC state (see Appendix [5] for the 
discussion) . 



V. DESCRIPTION ON ij-SECTIONS: UNITARY 
EVOLUTION 

In the case of the WDW quantization, we introduced 
in Subsec. IIII CI the notion of evolution by means of a 
family of observables which are related via unitary trans- 
formations. In this section we will analyze the possibil- 
ity of performing an analogous construction in the LQC 
model. First, in Subsec. IV A| we will establish the re- 
lation between the physical Hilbert space and an appro- 
priate space of "initial" data defined on a single slice 
Vi = const. That relation will be used in Subsec. IV Bl 
to construct a direct analog of the family (|3.14[) , which 
however fails to admit a unitary relation. In Subsec. IV CI 
certain modification of the construction will allow us to 
overcome this problem, although at the price of loosing a 
neat physical interpretation of the selected observables, 
which is recovered only in the large v\ limit. The modi- 
fied observables are finally used in Subsec. IVDI to extract 
physical predictions, which are presented in Subsec. IV El 



A. Time slices and associated Hilbert spaces 

Let us start with the general form of the wave function 
that represents the physical state, given by Eq. (|2.12p . In 
analogy with the procedure explained in Subsec. IIII CI 
we choose as the internal time the variable v\ and define 
the "initial data" functions on each slice v\ — const in 
the following way 



*«l(«2) V3) 



dw 2 dw3$. Ul (w 2 , u 3 )e 6 J 2 (v 2 )e £ J (v 3 ), 



(5.1) 
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where the spectral profiles ^ Vl {u>2^3) of $> Vl (i>2, U3) be- 
long to the slice Hilbert spaces 

H V1 :=L 2 (M 2 ,|^ 2 +^3||e^(«i)r 2 d^ 2 d W 3), (5.2) 

and are defined by the transformation P Vl : 7i? hy — > 7i Vl 



:= $(w 2 ,w 3 )e^ i (w 1 ). 



(5.3) 



On any slice «i = const, ej 1 (^i) provides just a func- 
tion of u>x which turns out to vanish in a set of zero 
measure [see Eq. (45) in Ref. 9] for the details]. There- 
fore, the map P Vl is unitary. This property, together 
with the fact that both e e J 2 (v2) and e^f (V3) form bases of 
their corresponding kinematical spaces Hf a , allows us to 
determine $ (0^2,^3) from $ Vl (v2,V3) (up to a zero mea- 
sure set). As a consequence, the projection on each v\- 
slice contains the same information as the entire physical 
solution. However, one cannot write the inner product 
of (|5.2p as an integral of <§> Vl («2, ^3) with well defined 
Lebesgue measure. Owing to the dependence of le^, 1 (vi)| 
in ll>2 and W3, the inner product of TL Vl is nonlocal when 
expressed in terms of V2 and V3. 

The unitary transformation P Vl allows us to define a 
map between initial data spaces. Each state on the phys- 
ical Hilbert space 7i? hy is associated, through P Vl , with 
a sequence of elements of the slice spaces Tt Vl . Each se- 
quence consists in the chain of "evolution steps" enumer- 
ated by v\ S . However, the corresponding evolution 
is not unitary, because under the identification of differ- 
ent slices, initial data belonging to one of the spaces TL Vl 
will in general not belong to the others, since (e^Vi)! 
depends on v% . In the next subsection we will try to pro- 
vide a more sophisticated notion of evolution free of this 
problem by building a set of observables analogous to the 
family (|3.14|) that we constructed for the WDW model. 



B. vi -observables 

Once we have introduced the Hilbert spaces (|5.2[) , and 
the transformations between them and 7i? hy , we can fol- 
low the construction of relational observables made in 
Subsec. IIII CI starting from the kinematical observables 
ln(t) a ) (a = 2,3), which also act as multiplication op- 
erators here. However, unlike the WDW eigenfunctions 
e a , i (ui), the eigenfunctions e^, 1 have a phase which 
is wi-independent separately on each of the subsemilat- 
tices and ^C^ i+2 , with a global phase shift of n/2 

between them (see Subsec. Ill Cp . Therefore the transfor- 
mation between TL V1 and Ti' analogous to (|3.12[) essen- 
tially removes all the information from the state. As a 
consequence, the observables constructed in this way do 
not carry physically interesting information. 

As an alternative, one may adopt a more naive ap- 
proach, which consists in considering the operators ln(u a ) 



just as multiplication operators acting on the elements of 
Ti vi in the w a -representation: 



[ln(v a )$ Vl ]{v2,v 3 ) = ln(v a )^ Vl (v2,v 3 ). 



(5.4) 



We can rewrite the action of these operators in terms of 
the variables w a and represent them as operators on the 
physical Hilbert space H? hy using Eq. (|5.3|) . In particu- 
lar, Mv 2 ) Vl acts on $ € S(R 2 ) C ?iS hy as follows: 

1 

[ln^m ^] (^2, 



e£\(fi) 



(5.5) 

whereas the action of In (O3) is the same with the subindex 
2 replaced with 3. This implies immediately that two op- 
erators at different times, e.g. \n(v a ) Vl : H^. hy — > H^ hy 

and \n(v a )vi : ^? y — > V > are related via the trans- 
formation 



w Phy 



$(cj 2 , w 3 ), 



so that 



(5.6) 



(5.7) 



Since the amplitude (i^i) | changes significantly 

both when u>x or v\ varies, the operators Q Vl . v , 1 are not 
unitary. Hence, the family of observables defined here 
fails to be unitarily related. 

In order to attain a notion of nontrivial unitary evolu- 
tion in vi, we propose in the next subsection a particu- 
lar construction which exploits the asymptotic properties 
of ej 1 (ui) and their relation with their WDW analogs 

C. ui-observables on components 



The success of the construction of Subsec. IIII Cl to pro- 
vide a nontrivial evolution picture for the WDW model 
rests in the form of the eigenfunctions (v 1 ) of the op- 
erator Gj, which are essentially rotating complex func- 
tions. In LQC, the analogous eigenfunctions e^Ui) os- 
cillate rather than rotate. Furthermore, on each of the 
subsemilattices and ( 4 lC+ +2 , these elements con- 

verge to a combination of incoming and outgoing WDW 
eigenfunctions, both contributing with equal amplitude. 
In this sense, each eigenfunction e% (vi) of ©1 can be in- 
terpreted as a standing wave, which contains both com- 
ponents moving forward and backward in time. This 
interpretation is supported by the studies of the classical 
effective dynamics of the system performed in Ref. [20| . 
where a bounce in the internal time v\ is observed. 
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These considerations suggest that, rather than try- 
ing to construct the analogs of the WDW observables 
]n(v a ) Vl , one should build instead two separate fami- 
lies ln(t) ){j: , each corresponding to one of the two com- 
mented components of the wave function. With respect 
to the procedure specified in Subsec. MI C| this can be 
viewed as a specific choice of two (instead of one) auxil- 
iary Hilbert spaces: TC' + and 7i'~ . 

In order to define the decomposition in a precise form, 
we first introduce the following transformation of the 
eigenfunctions of 1; defined in the distributional sense 
on each of the subsemilattices ( 4 ^C+ and ^C^ i+2 sepa- 
rately: 



T- 1 e[-s{b 1 -Ti/2)]Te 



(5.8) 



where s E {+, — }, e% S {e\,e\ + 2}, 9 is a Heaviside step 
function, b\ is the momentum conjugate to Vi [see Eqs. 
(|3.24[) and (|3.25[) ] , and T is a discrete Fourier transform 
analogous to the one defined for isotropic systems in Refs. 
Gill: 



Using these projections we can finally define two families 
of observables: ln(w )+ : Tt Phy — > H? hy and ln(« )~ : 
^Phy _^ -^Phy^ starting from the kinematical operators 
ln(t) a ). Their action on 7i? hy is analogous to Eq. (15. 5| . 



[ln(v a ) s Vl $](u) 2 ,u 3 ) 



d^ 2 ^ 



(«l) *(^2: W 3)- 



Within each particular family labeled by a and s, two 
observables evaluated at different times v\ and v* are 
related via the operators 



«C*(t*)*, 



[Q s vijVt $}(io 2 ,u; 3 ) = x;;,;; $(M2,a;a), 



(5.13) 



in the following way 



P7](&i) = E /(»o«i 



h e [0,tt]. 



(5.9) 

The introduction of the rescaling by v\~ 5 in this transfor- 
mation is needed to develop the analysis of the evolution 
in terms of b\ that we carry out in Sec. IVIl Therefore, 
we use the same Fourier transform here. 

The transformation (|5.8p essentially extracts in each 
subsemilattice the components (labeled by + and -) of 
efj that respectively converge, in the large v\ limit, to 



i, , i and e\, , i 



which move backward 



the WDW analogs e. 
and forward in time 

The functions e e J_^ sum up to the original eigenfunc- 
tions e^ 1 , therefore one can split any wave function $ Vl , 
defined in Eq. (|5.1|) . into rotating components simply 
by replacing the eigenfunctions e^, 1 in Eq. (|5.3p with e^ s 
[35j |. However, the Hilbert spaces D 3?^, which are 
the analogs of 7^ [in the sense of the definition (|5.2p ]. 
still have different inner products for different v±, since 
\e^"\ depends on v\. Thus, to "synchronize" the norms 
we include one more step in the splitting, namely the 
normalization of e^ s into pure phases, and introduce the 
corresponding auxiliary Hilbert spaces H.' s , analogs of 
(l3~TTp . The final splitting H^ hy -> H' s is thus defined 
(on the bases of the Hilbert spaces) as follows 



k>2 + ^3 



(5.10) 



The implementation of this splitting allows us to define 
the projection i?* of the physical states onto rotating 
components: 



R s ■ H P - hy 



H' s , H' s = L 2 (R 2 ,dLu 2 duj 3 ), 

[ii^ $] (wa, w 3 ) =xSi(w2»W3) 

:= l>(a;2,W3)e'^ 1 1 s (ui). 



H53HI 



ln(v a )l* =Q s Vuvt Hv a ) s Vi Ql 



(5.14) 



Since, by definition (|5.10p . | e'^ 8 (ui ) | = \u>2 2 for all 
v\ E Cj t , the operators Q s Vi v * are unitary and, therefore, 
within each family the considered observables are unitar- 
ily related. We can again extend the set formed by these 
families, adding the operators ©ali^ := ©a, to obtain a 
complete set of observables. 

Thus, the operators defined in Eq. (|5.12p provide a 
correct notion of unitary evolution. However, this comes 
at a price. Owing to the normalization (|5.10p . the ob- 
servables no longer have a precise physical interpretation. 
Such an interpretation can be recovered only asymptot- 
ically for large v\, where the rotating components e^ 1 / 
approach their WDW analogs (see Subsec. IIICp and 
[e^l converges to an w a -independent function [see Eq. 
(|3.4p ]. As a consequence, one can interpret the operators 
ln(?) a )* only approximately as evaluating ln(w a ) at the 
given value of v\ on the component that is moving for- 
ward (for negative sign) or backward (for positive sign) 
in time. The approximation improves as v\ increases; 
however, for v\ of the order of u>i or smaller (where the 
effective theory predicts a bounce in v\ ) all the precision 
is lost. An illustrative argument which shows the "unre- 
liability" of the interpretation of ln({i a )J in such regime 
is presented in Subsec. IV El where we discuss the ap- 
plication of the construction introduced here to analyze 
the dynamics of physical states which are semiclassical 
at late times. 



D. Numerical aspects of the analysis 

In this subsection we describe the numerical meth- 
ods used to analyze the dynamics of the model. The 
reader that is not interested in these numerical aspects 
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can safely skip this part and go to Subsec. IV El where 
the results are presented. 

We focus our discussion on physical states that are 
semiclassical at late times; more precisely, on Gaussian 
states (|4.14p peaked around large (oj^co^) [whose late 
time trajectory is determined by (/3 2 ,/3 3 ) via Eq. (|4.16|) ]. 
The wave functions corresponding to such states, given 
by the integral (|2.12p . are next evaluated applying the 
trapezoid method in the domain uj a £ [w*—5a a ,u*+5a a ]. 
$ (102,^3) has been probed within the uniform grid de- 
fined by the split of the domain into at least 2w* subin- 
tervals in each direction. 

To find the expectation values of ln(u )* we have used 
the expression of the elements x^ (^2,^3) of W s as func- 
tions of v a , 

X 8 Vl {v2,v 3 )= / dw 2 dw 3 x i Ul (w2,u;3)e^ 2 (-D2)e^(i;3). 

(5.15) 

Since eff (v a ) form orthonormal bases on their respec- 
tive kinematical spaces H.f , the physical inner prod- 
uct (|2.14[) on Ti ls takes a very simple form in the v a - 
representation 

(X S V1 \X' V \) = E X S V1 (V2^)X' V S M,V 3 ). (5.16) 

Furthermore, on ( u 2j V3) the observables ln^)^ act 
just as multiplication operators 

[Hva) s Vl X s Vl ](v 2 ,v 3 ) = ]n(v a )x s Vl (v2> v a)- ( 5 - 17 ) 
Hence, their expectation values on the state $ are 

(5.18) 

where Xv, ^ s related to $ via Eq. (|5.1ip . 

The dispersions corresponding to rn(u )* are given by 
the standard formula 

(A\n(v a )tf = (Qn^J 2 ) - (Hv a y vi }\ (5.19) 

where the expectation values of [ln(?) a )£ ] 2 are evaluated 
as we have explained for ($| ]n.(v a )„ 

In the numerical simulations, in order to calculate 
these expectation values we have first evaluated the com- 
ponents e^ s given in Eq. (|5.8[) with the use of the Fast 
Fourier Transform (FFT) algorithm, and then the nor- 
malized components e' £ J^ according to Eq. (|5.10p . Af- 
ter calculating the auxiliary profiles Xvii^,^), defined 
by Eq. (|5.11[) . we have evaluated their wave functions 
Xvx (^2,^3) via direct integration of Eq. (|5.15p , in a simi- 
lar way as explained for the evaluation of Next, we 
have found the expectation values computing the sum- 
mation in Eq. (|5.18p . The domains of v a have been 
restricted to £+ n [0,4a;*]. The values of v\ have been 
chosen from the set n [0, 2 • 10 5 ]. Finally, we have 
performed simulations for various values of e and for lu* 
ranging from 2.5 • 10 2 to 10 3 . 



E. Results and discussion 

In Fig. [5] we show two examples of the results of our 
numerical study. In all cases, the analysis reveals that 

• Inasmuch as the expectation values (|5.18p are con- 
cerned, the states remain sharply peaked for all val- 
ues of v\ . This applies to both components moving 
forward and backward in time. 

• For v\ 3> u>*, the expectation values follow the clas- 
sical trajectories, whereas, when the universe ap- 
proaches the classical singularity, the discrete ge- 
ometry effects induce repulsive forces which cause 
bounces in both V2 and 113 at the values predicted 
by the effective dynamics. 

These results prove the robustness of the big bounce 
scenario of LQC, extending its validity to another system: 
the vacuum Bianchi I cosmologies. They confirm as well 
the ability of the effective dynamics to predict correct 
results with errors much lower than the spread of the 
wave function. 

Nonetheless, there is one aspect of the results which re- 
quires special comments. As one can see in Fig. [5J while 
for large v\ the expectation values of the observables 
ln^o)* follow the classical trajectories (with decreasing 
precision for decreasing V\), for v\ < 0.2 wi(w*, u;*) the 
expectation values of both families "freeze" at the same 
trajectory of constant V2 and ^3. However this behav- 
ior does not correspond to that of ln(ii a ) at given v\ in 
any physical way, since \n.(v a )% no longer approximates 
it (for each a) in that region. 

To illustrate this fact, let us consider a slight modi- 
fication of the proposed scheme for the construction of 
observables. So far, we have used v\ as an emergent 
time; however, owing to the symmetry of the system, 
any of the variables «j can play this role. Therefore, for a 
given physical state which is semiclassical at late times, 
we can consider, e.g., the two following constructions: 
\n.(v a ) s Vi and ln(« a ')J , where a 1 = 1,3, and compare the 
corresponding expectation values. Such comparison is 
presented in Fig. [31 Outside the regions where the effec- 
tive dynamics predicts bounces, the trajectories formed 
by the two considered families agree. However, once the 
value of the emergent time corresponding to each of the 
families drops below that of the point where the effective 
bounce occurs, we observe significant discrepancies in the 
resulting trajectories. Indeed, while one family shows a 
bounce, the other predicts the freezing of the trajectory, 
and vice versa. This implies that in such regions the tra- 
jectories cannot be interpreted as corresponding to the 
value of ln(i>i) at given emergent time in any (even ap- 
proximate) sense. The applicability of the observables 
associated with Ui-slices and with Vi as an emergent time 
is thus limited from a physical viewpoint. 

In order to obtain a more reliable description, with a 
valid physical interpretation in all the interesting regions, 
we have to use a different construction of observables. 
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FIG. 2: For Gaussian states [with profile (|4.14p ], the expectation values and dispersions of the observables lnfOa)^ and 
lnfOo)^ [corresponding to epochs when the wave packet is moving forward (green errorbars) and backward (red errorbars) in 
time, respectively] are compared with classical (pink lines) and effective (blue dotted line) trajectories. The states are peaked 
at u>2 = <^3 = 10 3 with the relative dispersions AO2/O2 = AOa/03 = 0.05. The phases are, respectively, /3 2 = /3 3 = for 
(a), and /3 2 = /3 3 = 0.1 for (b). The expectation values follow the effective trajectory till the point of bounce in v\, where 
the dynamics "freezes". On the other hand away from the bounce the quantum trajectory approaches the classical ones. In 
particular, in (a) the trajectories before and after the bounce coincide on the plane Vi — «2- 
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FIG. 3: Comparison of the expectation values of the two fam- 
ilies of observables ln^)^ (blue) and ln(«i)+ 2 (red), corre- 
sponding to two different choices of emergent time, calculated 
on the state considered in Fig. [2^. 



Such construction will be presented in the next section, 
where as emergent time, instead of v±, we choose its con- 
jugate momentum b\. 



properties of the eigenfunctions (see Subsec. Ill C[) shows 
that v\ is not monotonous in the proper time (see also 
Appendix [B]) . This forced us to define two families of 
operators corresponding to epochs when the universe is 
respectively contracting and expanding in v\. Such a 
splitting can be avoided if we identify a phase space vari- 
able which changes monotonously with time. Here again 
the effective dynamics comes to the aid, showing that any 
canonically conjugate momentum hi of Vi possesses the 
desired property, so that it may be a promising candidate 
for an emergent time in the genuine quantum theory. 

In this section we explore this possibility by replacing 
the configuration variable v± with its momentum b\ and 
representing H? hy by Hilbert spaces of data on 61 = const 
slices (in analogy with the v± -slices discussed in Subsec. 
IV Aj) . This step is described in Subsec. IVI Al We then 
introduce in Subsec. IVI Bl the observables corresponding 
to ln(i>i) at a given b±. The constructed family of oper- 
ators is used in Subsec. IVI CI to carry out a numerical 
analysis of physical states, whose results are presented in 
Subsec. IvTdI 



A. 6-representation 



VI. UNITARY EVOLUTION ON 6-SECTIONS 

One of the main reasons why the analog of the WDW 
observables defined on v\ -slices fail to provide a unitary 
evolution in the loop quantization is the fact that the 
variable chosen as emergent time does not possess all the 
ideal properties of this kind of time. Indeed, an analy- 
sis of the dynamics at the effective level [2(| and of the 



Let us start by choosing at the classical level the mo- 
mentum bi conjugate to v\ which is defined through Eqs. 
(|3.24p and (|3.25p . In our quantization, the transforma- 
tion between the "position" and "momentum" represen- 
tations is given by the formula 
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As commented in Subsec. IV C[ we have introduced an 
additional rescaling by | x / 2 in comparison with Eq. 
(|3.26p . The reason for that is twofold: the rescaling 
guarantees a proper convergence in order to perform the 
transformation numerically and, in addition, ensures a 
more convenient behavior of the transformed eigenfunc- 
tions (see Subsec. IVI Cp . Under this transformation, the 
kincmatical operators which enter the Hamiltonian con- 
straint and the observables change as follows 

Vi — > 2id bi , 

H 1/2 (#2* -#-2*) \vi\- 1/2 -> 2ism(bi), (6.2) 

where 6, acts in the new representation by multiplication. 

Since, in order to construct the physical state, we re- 
strict the analysis to one superselection sector, with func- 
tions supported on , the same restriction applies to the 
transformation (|6.ip . In this case, the domain of bi can 
be taken as the unit circle bi G S 1 . Moreover, to further 
adapt the transformation to our specific system, let us 
note the following: 

(i) The restrictions of the cigcnfunctions of 8^ (sup- 
ported on £+) to the subsemilattices and 

£*t-+2 are eigenf unctions of 0f , which differs from 
the evolution operator of the isotropic system just 
by a compact operator (see Ref. [9j and Subsec. 

BE). 

(ii) Since, e.g., any v\ = const section contains all the 
information about a physical state, in particular, 
the restriction of a physical state to one of the sub- 
semilattices for v\ uniquely determines the other. 



where e* 1 (&i) is 



sum- > ^*»> 



bi e [0,tt) 



[^ae&K&i-Tr), 6iG[7r,27r), 



(6.5) 



As discussed above, we restrict our considerations, for 
example, to the domain of b\ which corresponds to the 
first of these cases: b\ <E [0,7r). Within this domain, one 
can now apply the analog of the constructions presented 
in Sec. IIIII and Sec. [V] In doing so, we introduce the 
counterpart of the spaces H Vl [given in Eq. (|5.2[) ]. which 
are defined on surfaces of constant b\, 

H bl ~L 2 (R 2 , \uj 2 +LJ 3 \\el\(h)\- 2 dLJ 2 duj 3 ), (6.6) 

and the transformation of the physical states into ele- 
ments in these spaces 



(6.7a) 
(6.7b) 



Each function <& bl at any given slice b\ = const deter- 
mines the physical state uniquely, thus it can be under- 
stood as the initial data at that slice. Therefore, one 
may develop a notion of evolution in b\ using the cor- 
responding map between initial data spaces. It is worth 
noticing, nonetheless, that, unlike in sLQC (and with the 
chosen factor ordering) , the Hamiltonian constraint (|2.5|) 
in a ^-representation cannot be expressed as a differen- 
tial equation of finite order. 



(iii) In the isotropic model and within the context of 
solvable LQC (sLQC) [12| a similar transformation 
allows one to select the domain of the corresponding 
momentum b as [0, 7r). Furthermore, the sLQC ana- 
log of 2 is proportional to [sin(6)9b] 2 . This in turn 
allows one to write the constraint as a Klein-Gordon 
equation. 

This last point indicates that it is convenient to apply 
the transformation (|6.1[) to functions supported on ^Cf. 
and ' 4 '£+ +2 independently. The domain of bi consists 
then of two copies of the circle with radius one half (i.e., 
the periodicity of bi is 7r in each copy). As a consequence, 
one can define the transformation T such that 

[^](6i)= i>(vi)vi-h-i Vibi , (6.3) 

Vi £(*)£+ 

and, analogously, the transformation T 2 by choosing the 
summation over ^C^. +2 instead of 

If we transform the cigcnfunctions of Oi, we can define 
a "mixed" representation of physical states 

®(h,v 2 ,v 3 )= du> 2 dw3$(u;2, u 3 )e s J (bi) 

Jr2 (6.4) 

x e £ J 2 {v 2 )e £ J 3 (v 3 ), 



B. Observables associated with ^-slices 

With the spaces of initial data at constant b\ at our 
disposal, we can now define the observables ln^)^. For 
i = 2, 3 the construction is completely parallel to that 
presented in Subsec. IV BI however, it suffers from a sim- 
ilar problem of lack of a unitary relation between the 
operators in each family of observables. Here, neverthe- 
less, the problem is much less severe. Indeed, the func- 
tions 13^(61) are (except possibly for small u)\) rotating 
functions that never vanish (see Subsec. IVI C[) . This 
property has a qualitative analytical explanation. Ac- 
tually, the operator T\{Q\) is, up to a compact correc- 
tion, equal to — [127T7G sin(6i)9h 1 ] 2 , whose eigenspaces 
are spanned by basis elements of the form jV"e ± * feB ' 6l ) ) 
with x(bi) := ln[tan(&i/2)] and where N = N(k) is a 
normalization factor. Furthermore, once we restrict the 
space to functions supported on v\ > (as happens to 
be the case in our model), the eigenfunctions reduce just 
to purely rotating functions [with vanishing contribution 
of one of the two phases ±ikx(bi)]. These eigenfunctions 
will differ from those of the actual operator just 
by small corrections whose (kinematical) norm decreases 
with uj\. 
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As a consequence, we do not need to split the eigen- 
functions into rotating components in order to form a 
family of unitarily related observables. The following 
parallel of the transformation (|5.10[) : 

i e £1 (b) 

el\(b) - e%(b) = \u 2 +^jg±± (6.8) 

is sufficient to build analogs Rb 1 of the operators R* , 

R bl : -> H' bl := L 2 (R 2 , dw 2 dw 3 ) , (6.9a) 

[R bl $}{u 2 ,uj 3 ) = x bl (w 2 ,^3) (6.9b) 
:=$(w 2j w 3 )e , S L I (6i). 

The observables ln(i)i)b 1 are defined as follows. For 
i = 2, 3 we simply act on 

Xb 1 (v2,v 3 ) = / du 2 dui 3 x bl (uj2,uj 3 )e £ J 2 (v2)e £ J(v 3 ). 

JR 2 

(6.10) 

with the kinematical operator ln(uj). Thus, the action of 
ln(w2)b! on the elements of 7iS hy turns out to be 

x | d^(e-Jm(« 2 )e-) K + e^^^^&O^C^,^), 

and similarly for i — 3. These families are unitarily re- 
lated by the very same reason explained for ln(v a ) Vl in 
Subsec. IV C| and can be extended equally to a complete 
set of observables. 

As an aside, let us comment that in the sLQC approach 
the analog of the basis eigenfunctions e^, 1 (&i) are pure 
phases, and then the transformation (|6.8[) reduces just 
to a unitary rescaling by \u>2 + uj 3 \^ . 

The construction of hi(vi)b 1 at a heuristic level is based 
again on the action of the kinematical operator on some 
suitable modification of the initial data at b\. To give a 
precise definition, we first need to express the operator 
ln('0i)b 1 in the 6i-representation. 

Note that the function ln(ui) can be expanded around 
a point v° > as 

oo 

ln(V) = mK)-£-— T-M-Vt)" (6.12) 
i n[v? ) n 

n—l v 1 ' 

Promoting the terms u" to operators and applying rela- 
tion l|6.2p . we can represent the operator ln({Si) as 

\xx{v x ) = ln«) - "T^wK ~ 2zd bi r, (6.13) 

n— 1 v 1 ' 

defined on elements \ bl ( v 2, v 3 ) of the kinematical Hilbcrt 
space. Changing to the w a -representation via Eq. (|6.10D . 
the action on % bi (w2,W3) can be expressed as 

NtfOXfeK^.ws) = H"2^3)Mn)e%]{bi), (6.14) 



and thus the final form of these observables reads: 
lnfa)* : y -> W £ ? hy , 
pn(«iK*](wa.«3) = [ln( !4'Sr i) |> (^^3),(6.15) 

with ln(-Ci) given in Eq. (|6.13[) . 

One can check by inspection that these observables are 
not related by a unitary operator Qb u b* '■ T i^ hy — > W? hy , 
analogous to the operators defined in Eq. (|5.13|) . There- 
fore, the observables ln({>i)h 1 can play only an auxiliary 
role with respect to the families of unitarily related ob- 
servables ln^a)^. In the next section we employ these 
families of observables to analyze the dynamics of phys- 
ical states which are semiclassical at late times. 



C. Numerical analysis 

Let us briefly discuss the numerical methods used in 
our calculations. This is a technical subsection that can 
be safely skipped if the reader is not interested in the 
numerics. 

As in Sec. [V] we concentrate our study on the Gaus- 
sian states (|4.14[) . In most cases, the calculations are 
exactly the same as those of Subsec. IVDI for the observ- 
ables associated with wi-slices. In particular, to compute 
the expectation values and dispersions of ln(w a ) we have 
used the analogs of Eqs. (|5.15|1 - (|5.19[) . replacing x s Vl with 
X bl ■ The eigenfunctions e'^ 1 . , needed to determine x bl > 
have been computed via Eq. (|6.5[) by means of a FFT 
algorithm. In comparison with Subsec. IV Dl the only rel- 
evant difference appears in the calculation of (ln({ii)f )1 ) 
and (Aln({ii)h 1 ), which is performed as follows. 

In a first step, we calculate the action of ln(vi)b 1 on $ 
using Eq. (|6.15jl . To find the values of lB.(vi)e'^ , we act 
with ln(£>i) in the wi-representation, and transform the 
result to the ^-representation via Eq. (|6.1[) by applying 
a FFT. 

At this stage, however, it is necessary to comment on 
a technical subtlety. In order to calculate the desired re- 
sult exactly as specified in the previous subsection, one 
has to normalize the eigenfunction according to Eq. (|6.8p 
before acting on it with hi(vi). Since the normalization 
must be done in the 61-representation and the action of 
the operator is known instead in the wi-representation, 
one would have to carry out a sequence of operations, 
namely, a Fourier transform, a normalization, and an in- 
verse Fourier transform. Unfortunately, the FFT algo- 
rithm used for this assumes that both the function that 
is transformed and the result are supported on points 
which are uniformly distributed in the circle. Since, in 
our case, the support of the eigenfunction is an entire 
subsemilattice, ^Cf , one should restrict it to a set of 

the form ^Cf R [0,i> max ] (where v max has some large 
value), thus removing the tail, which decays (up to loga- 
rithmic factors) as v^ 1 . The restriction would give rise to 
numerical errors which would manifest themselves near 
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bx = and b\ = tt. In turn, this would lead to errors 
in normalization, which would be amplified through the 
commented sequence of operations (FFT — » Normaliza- 
tion — > FFT), making the results in the large volume 
regime unreliable. 

In order to avoid this problem, we implement a slightly 
different algorithm. Instead of normalizing e%J- , we first 
act with ln('Ci) on the nonnormalized eigenfunctions, car- 
rying out the normalization afterwards. That is, 



(6.16) 



Such a method introduces an error owing to the differ- 
ence between the normalized and nonnormalized func- 
tions. However, one can argue that the error can be ne- 
glected for semiclassical states peaked around large u*, 
using that e'^, 1 can be approximated by its sLQC 
counterpart in that regime: 



~ |W2+W3|* [e Mll(il, +OK- 2 ) 
x(b) := In 



, b 

tan - 
2 



(6.17) 



Once the functions are known, we evaluate the 
profiles 



[M^i)XhJ( w 2, v 3 ) = / dw 2 dw 3 ^(w2,w 3 )/^(6i) 

JR 2 

x ^ 2 > 2 )e5> 3 ), 

integrating in practice the righthand side over the inter- 
vals u a £ [uj* — 5er*, uj* + 5er*] via the trapezoid method. 

The computed profiles are finally used to calculate the 
expectation values 

($|ln(« 1 ) fel $) = || Xbl ||- 2 (EH]) 

X ^2Xb 1 ( V ^^ V 3)[ ln (vi)Xb 1 ]( v 2,V 3 ), 



D. Results 

A representative example of the numerical analysis of 
the Gaussian states is presented in Fig. |U The general 
conclusions of this analysis can be summarized as follows: 

• The states that are sharply peaked at some initial 
b\ remain so during the whole evolution. 

• For small bi, the expectation values of the three 
families of observables ln({)j){, 1 follow classical tra- 
jectories corresponding to a universe contracting in 
the three directions. As b\ increases and the values 
of Vi decrease becoming of the order of w*, we ob- 
serve a deviation from the classical dynamics which 
results in independent bounces in all the three di- 
rections. After each bounce (in direction i), the 
value of Vi starts growing, and its dynamics quickly 
approaches that of a classical expanding universe. 

• Through the entire evolution, the genuine quantum 
trajectories agree to a good precision (less than 10% 
of dispersion) with those corresponding to the clas- 
sical effective dynamics presented in Appendix IBl 

The results listed above provide an additional (inde- 
pendent) confirmation of the conclusions presented in 
Subsec. IV El However, the current analysis constitutes a 
significant improvement with respect to that of Sec. [V] 
Namely, while it is still only in an asymptotic sense that 
the observables used here possess a well defined physi- 
cal interpretation, the accuracy of this interpretation re- 
mains fairly good during the entire evolution. Unlike the 
situation found for hx(v a ) Vl , where the observables com- 
pletely lose a meaningful physical interpretation in some 
epochs of the evolution, the level of confidence depends 
now just on the peak values u>*. Furthermore, this level 
improves as one considers more and more macroscopic 
universes (larger to*). 



VII. CONCLUSIONS AND DISCUSSION 



where owing to technical limitations the summation was 
restricted to C 2 := (£+ n [0,44]) x (£+ n [0,44]). 

The dispersions of these observables are found using 
the standard relations [similar to Eq. (|5 . 19[) ] . The ex- 
pectation values ($| kx(vi)% $>) are computed with an 
algorithm that is fully analogous to the one used for 
($|ln(t>i)6 1 *). 

In our numerical simulations, the number of points se- 
lected for the FFT of the eigenfunctions is equal to 2 17 , 
distributed uniformly in the set b\ € [0, 7r). The simu- 
lations have been performed for the same range of pa- 
rameters as in Subsec. IV Dl that is, for ui* ranging from 
2.5 • 10 2 to 10 3 , and for relative dispersions in u a be- 
tween 0.05 and 0.1. The results of these simulations are 
discussed in the next subsection. 



In this work we have investigated the concept of evo- 
lution in LQC. Our main goal has been the construction 
of a well defined unitary evolution framework with no 
reference to a matter field as a clock. We have achieved 
this by building certain families of unitarily related par- 
tial observables, parameterized by geometry degrees of 
freedom. As a test bed for our study, we have chosen 
the specific example of a vacuum Bianchi I toroidal (T 3 ) 
universe. 

Within this model, we have proposed two different con- 
structions of observables, built out of the kinematical op- 
erators ln(v a ) with a = 2,3, and parameterized respec- 
tively by: (i) the affine parameter v% corresponding to 
one of the triad coefficients p\ , which plays the role of a 
configuration variable, (ii) its conjugate momentum b±. 
In this way, the role of internal time has been assigned 
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FIG. 4: Expectation values and dispersions of the observables ln^j)^ on the Gaussian states [with profile (|4.14|l ] presented in 
the V2 — Vi plane and compared with classical (green lines) and effective (blue dotted line) trajectories. The states are peaked 
at u>2 — ui* = 10 3 with relative dispersions AO2/02 = AOs/03 = 0.05. The phases are, respectively, (3 2 = /3 3 = for (a) 
and ft 2 — A — 0.1 for (6). The expectation values follow the effective trajectory through all the evolution. In particular, the 
trajectories before and after the bounce coincide in (a). 



to geometry variables, making its definition independent 
of the matter content. 

In order to develop and test our constructions in an 
analytically solvable setting, as well as to identify the im- 
print of the loop quantization in the dynamics, we have 
first implemented them in the quantum version of the 
system obtained via the standard methods of geometro- 
dynamics, i.e., in the WDW analog of our LQC model. 
In both of the LQC and WDW theories, the physical 
Hilbert space of our model is L 2 (W 2 , \u>2 + o->3\duj2<iui3) . 
The difference between them is captured in the form of 
the observables. 

In the WDW theory the basis elements of the kinemat- 
ical Hilbert space are, essentially, rotating functions, and 
owing to this property the unitarity in the evolution is 
easily achieved both for cases (i) and (ii). Furthermore, 
all the constructed observables have a neat physical inter- 
pretation as the evaluation of hi(v a ) at a given moment 
of vi or 61, respectively. 

The constructions described above have been used to 
investigate the dynamics of states which are semiclassi- 
cal at late times, with a sharply peaked Gaussian profile. 
We have shown analytically that any such state follows 
the classical trajectory right to the classical singularities, 
which arc thus not resolved. Actually, this lack of singu- 
larity resolution has been shown to be a general property 
of the physical states, not only of the Gaussian ones. 

In the case of an LQC model, the situation is more 
complicated. While in case (ii) we have succeeded in 
applying the construction tested in the WDW model al- 
most straightforwardly, case (i) has posed a considerable 
challenge. Both a direct application or a straightforward 
generalization of the scheme introduced for WDW have 
resulted in operators which either fail to provide a uni- 
tarily related family of observables, or do not encode any 



physically interesting information about the system. 

To overcome this problem, we have used our knowl- 
edge of the WDW limit of an LQC state, observing that 
any limit of this kind contains two types of components, 
which respectively move forward and backward in the 
evolution parameter. This feature has motivated the in- 
troduction of a similar splitting in rotating components 
of the exact LQC wave functions, which has been de- 
fined in a precise manner exploiting the properties of the 
Fourier transform. The scheme for the construction of 
observables defined in the WDW case has been applied 
to each of the components separately and has allowed us 
to attain the required unitarity in the evolution. 

In contrast to the WDW case, where the constructed 
observables admitted a neat physical interpretation, in 
LQC, owing to the modifications in the construction nec- 
essary to ensure unitarity, the physical interpretation of 
the operators is clear only in approximate sense. In par- 
ticular, for case (i) the interpretation of the operators 
as corresponding to the value of \n(v a ) at a given v\ is 
valid only in the large v± limit, and accurate only for 
part of the evolution. For Gaussian states peaked around 
some uj*, the interpretation ceases to be acceptable when 
v\ < 0.2 |cl>i(cl>2, WjOI, a feature that we have seen by ex- 
changing the spatial indexes of the kinematical operators 
and of the internal time v\. Nonetheless in case (ii), while 
the interpretation is still approximate, its accuracy does 
not change during evolution (i.e., it does not depend on 
the value of 61) and is pretty good for the physically in- 
teresting states (peaked at large uj*). 

Both constructions, (i) and (ii), have been employed to 
analyze the dynamics of Gaussian states sharply peaked 
around large uj*. In case (ii), and in order to bring the 
analysis to a common setting (v\-v a plane) to facilitate 
the comparison with (i), an additional family of observ- 
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ables -corresponding to ln(«i) at given 61- has been con- 
structed. However, these observables are not unitarily 
related; therefore they play only an auxiliary role in the 
discussion. The numerical analysis reveals the following 
picture (valid both for case (i) and (ii), unless otherwise 
stated): 

• The considered Gaussian states remain sharply 
peaked (as far as the relative dispersions of the con- 
structed observables are concerned) through all the 
evolution. 

• For large Vi/u>i(u>2, W3), the expectation values of 
the observables follow classical trajectories, while 
for other values of this ratio we have observed de- 
viations which ultimately lead to bounces in v a . 

• Since the physical interpretation of the observables 
in case (i) is not reliable for small v\, because 
Vi/u>i(u>2, k/3) becomes small, we have also per- 
formed the analysis in terms of the family of ob- 
servables of case (ii). This analysis has shown as 
well the bounce in v\. 

• The expectation values of the introduced observ- 
ables (calculated in the domains where their phys- 
ical interpretation is valid) agree with the trajec- 
tories predicted by the classical effective dynamics 
presented in Appendix [Bl 

As a consequence, in LQC the singularity, already re- 
solved at a kinematical level, is also resolved dynami- 
cally. Furthermore, the above observations imply that, 
similarly to what occurs in the isotropic system, large 
semiclassical universes expanding in terms of the coef- 
ficients pi are connected to also large and semiclassical 
contracting ones via a sequence of bounces which are in- 
dependent for each direction. 

Let us point out that the complication of the numer- 
ical problems that one has to solve in the analysis of 
the vacuum Bianchi I model has not allowed us to probe 
states as semiclassical as the ones studied in the isotropic 
case (considered for example in Ref. Q), nor check re- 
ally the long term evolution of the system. Besides, the 
numerical analysis has been restricted just to Gaussian 
states. This analysis alone is not sufficient to conclude 
that all states that are semiclassical before the bounces 
stay so after them. However, the semianalytical study 
of the WDW limit that we have performed for case (i) 
has shown that, whenever the state has finite expecta- 
tion values and dispersions for the set of operators (|4.27|) 
[and finite expectation values for (|4. 28[) ] , the relative dis- 
persions of ]xi(v a ) Vl in both branches of the evolution are 
asymptotically equal. 

In summary, taking the model of the vacuum Bianchi I 
universe as an example, we have constructed a notion of 
unitary evolution in a well defined way and without any 
reference to a matter field as a clock. We have achieved 
this goal by building several families of unitarily related 
observables, which in turn have allowed us to thoroughly 



analyze the dynamics of the considered model. The re- 
sults of this analysis confirm at a genuinely quantum 
level the robustness of the bounce picture, as well as the 
preservation of the semiclassicality across the bounces. 

To conclude, we would like to comment that the ap- 
proach presented here pioneers the development of a new 
methodology which can be employed whenever one has 
to analyze the evolution in absence of degrees of freedom 
that are quantized in a standard way. Therefore, it can 
be applied in a large variety of polymerically quantized 
systems. In particular, it is applicable to the quantiza- 
tion carried out in Ref. [10] where, for the analogs of the 
constructions (i) and (ii), we can use the total volume 
v defined in that reference and its conjugate momentum 
b, respectively. Nonetheless, the prescription for the so- 
called improved dynamics chosen for our study here offers 
us a precise control over the introduced constructions, 
inasmuch as the mathematical structure of the system 
has already been carefully studied in the literature and 
is known in great detail. Furthermore, it allows to test 
explicitly the reliability of the results provided by the 
observables. In particular, in the case of the operators 
parameterized by v\, we have been able to recognize a 
loss of predictability (at least in terms of a conventional 
physical interpretation of the observables) by interchang- 
ing the roles of internal time and dynamical variable be- 
tween v\ and v a (a — 2 or 3), a possibility which is not 
available in the prescription followed in Ref. (loj . 
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APPENDIX A: DISPERSION IN THE 
WHEELER-DEWITT LIMIT: ENSEMBLE OF 
SECTORS AND FURTHER CONSIDERATIONS 

In this appendix, we analyze the behavior of the disper- 
sions in the WDW regime taking into account the eight 
sectors in which the LQC wave functions have a well de- 
fined limit, rather than considering each of these sectors 
independently. We want to define the ensemble of those 
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sectors in such a way that the properties of the expecta- 
tion values and dispersions of the LQC observables are 
reflected by features of their analogs on that ensemble. 
To achieve this goal, we note the following properties of 
the operators lxi(v a ) Vl : 

(a) Since each operator is defined for a definite value of 
V\, which belongs just to one of the subscmilatticcs, 
( 4 ) C €l or ( 4 )£ £l+2 , we can consider the states sup- 
ported on them separately. As a consequence, we can 
safely consider two independent limits corresponding 
to the splitting of the support in v\ . 

(b) For given v\, the action of the operator \n(v a ) vi can 
be represented as the action of the sum of the respec- 
tive restrictions to each of the four sectors obtained 
with the splitting into subsemilattices in the direc- 
tions a = 2,3. 

Applying the calculation method explained in Subsec. 
IV Dl to the expectation values and dispersions, one sees 
immediately that 

(i) the expectation values of hi(v a ) Vl are arithmetic av- 
erages of the expectation values of the restrictions, 
and 

(ii) the dispersions are bounded by an arithmetic aver- 
age of the dispersions of each of the four contribut- 
ing sectors plus terms of the form \J (-D^), where 

the operator D m measures the difference between 
the total expectation value and that of the consid- 
ered sector, labeled by m = 1, ...4. 

On the other hand, in Subsec. IIV CI we proved that, 
for each of the sectors, the relative dispersion of ln(v a ) Vl 
(a = 2, 3) is the same in the two components moving for- 
ward/backward in time, respectively, provided that some 
mild conditions (e.g.) on the forward- moving component 
are satisfied. In order to extend the result to the en- 
semble of sectors, one has to cope with the possible dif- 
ferences between the expectation values of hi(v a ) Vl men- 
tioned above in point (ii). 

To verify the boundedness of these differences we ob- 
serve that the relation between the restrictions of the 
eigenfunctions to the subsemilattices corresponding to 
the different sectors under consideration -which is en- 
coded in Eq. (40) of Ref. 0]-, together with the relation 
between the normalization of these restrictions and that 
of their WDW limits, imply that the coefficients A a and 
W a are equal for all the studied sectors [see Eqs. (|3.19ap 
and (|3.22a[) ]. From the invariance of these coefficients 
under the transformation $_ — > $ + and the considera- 
tions made in Subsec. IIV CI it follows then that the an- 
alyzed differences have indeed a well defined, finite limit 
as V\ — ► oo. 

We also note that, as in the case of the discussion of 
Subsec. IIV CI the above arguments can be repeated con- 
sidering the transformation $ + — > opposite to the 
one that we have studied here. 



Finally, it is worth commenting that the requirement 
that the state have finite expectation values and disper- 

(2) 

sions with respect to the operators Q, a , although provid- 
ing a restriction on the space of states, can be considered 
quite reasonable from a physical viewpoint. Actually, 
since the choice of v\ as an emergent time (and hence 
the role of lo\ as a "time frequency") is arbitrary, one 
can exchange the roles of u)\ and u)% and, for a given 
state represented by $(^2,^3), consider its transforma- 
tion into a different one 

$ _j. $' : $'( W2 , oj 3 ) = $(^1(^2,^3)^3) (Al) 

On the other hand the transformation 

$(lu 2 ,lu 3 ) ■-> ^^^ l^^,^),^) (A2) 

corresponds just with swapping the coordinate v\ for v 2 . 
Similar transformations can be defined also for the in- 
terchange V\ <-> v 3 . Combining both of them one can 

rewrite the square of £l a as the square of ln^l act- 
ing on the state represented by Therefore, since that 
square is the main component of the dispersion, one can 
relate (A£l a ) with the dispersion of the operator In |w a | 
on $' at least heuristically. Similar arguments can be 
applied as well to £{, . 

APPENDIX B: EFFECTIVE CLASSICAL 
DYNAMICS 

The quantum dynamics resulting from the LQC ap- 
proach turns out to be surprisingly well reproduced by 
certain classical effective dynamics, constructed by re- 
placing the holonomies and fluxes in the quantum Hamil- 
tonian constraint by their expectation values [13, EH- 
This in principle naive construction has been proven to 
accurately mimic the genuine quantum evolution (with 
discrepancies much smaller than the quantum disper- 
sions) in several situations (see e.g. d, Q, HH). We 
sketch here its derivation for the case of the vacuum 
Bianchi I model. 

We start with the quantum constraint (|2.5p . applying 
to it the replacements explained above. The resulting 
effective Hamiltonian takes the form 

H cS = - \ [QfQf + Of Of + ©f Of ] , (Bla) 
Qf :=6w-fGviSm(bi), (Bib) 

where the variable hi (i = 1, 2, 3) is the momentum con- 
jugate to vi [see Eqs. (13.24[) and (]3 . 25[) ] . This momentum 
is compactified as the unit circle. 

Given the Hamiltonian (|B1[) . one can easily derive the 
Hamilton- Jacobi equations for d T Vi and <9 T 6 i; where r 
is a time parameter associated with the Hamiltonian. 
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Defining the constants of motion 9° ff = 8njGlCi and 
JC = K-i + JC 2 + /C 3 (as in Ref. [13]) we have 



(d T Vi) 2 = 9(8ttG) 4 (/C-A: 4 ) 2 

d T bi = 3(8ttG) 2 (/C j - £) sin(6,). 



4JQ 
3 



(B2a) 
(B2b) 



Equation (IB2b[) immediately implies that, within the 
intervals bi € [0, ir) and hi S [7r, 27r), which in the LQC 



quantization correspond to the respective subsemilattices 
( 4 )£+ and ( 4 lC+ +2 (see Sec. I*VT|) . each of the variables 6, 
provides a good candidate for an internal time. On the 
other hand, the general solution to Eq. (|B2a[) reads 



Vi(r) = -/C iC osh[3(8^G) 2 (/C - Ki)(r - r Q )], (B3) 

where t d is a constant representing the time of the 
bounce. Therefore, Vi is not monotonous in time. 
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